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Abstract 

Eynard's formulation of Hermitean 1-matrix models in terms of 
intrinsic quantities of an associated hyperelliptic Riemann surface is 
rephrased as a Lagrangean field theory of a scalar particle propa- 
gating on the hyperelliptic surface with multiple self-interactions and 
particle-source- interactions. Both types of interaction take place at 
the branch points of the hyperelliptic surface. 



1 Introduction 



The conceptual and practical understanding of some random matrix models 
has been advanced significantly by the seminal work of Eynard [4J. This 
author introduced an elegant new approach to the topological expansion of 
Hermitean 1-matrix models based on the intrinsic notions of a hyperellip- 
tic Riemann surface which is — as spectral surface — attached to the matrix 
model. The new approach has been further elucidated in [5] and applied to 
2-matrix models in [6] and [7j. It has moreover been shown in [8] that one 
may abstract from the framework of matrix models to the construction of 
new invariants attached to more general Riemann surfaces (mimicking the 
topological expansion of matrix models). This brings into focus jH] new sub- 
jects as Kontsevich's matrix model for intersection numbers [9], topological 
string theory [TQ] , [TT], [T2] and [13] . as well as further applications to math- 
ematical subjects, as the recursive determination of Weil-Petersson volumes 
[H], [T3] and a conjectural matrix model representation of Hurwitz numbers 

The purpose of this note is to show that Eynard's elegant construction can be 
deciphered into a non-elegant effective Lagrangean formalism. We will con- 
centrate here on Hermitean 1-matrix models — but we hope that our methods 
will also apply in some of the instances mentioned above. The Lagrange to 
be constructed will be that of a scalar field propagating on a hyperelliptic 
Riemann surface (i.e. the Riemann surface given as the spectral curve asso- 
ciated with the matrix model). The corresponding 'vacuum graphs' of the 
Gell-Mann-Low expansion and the graphs for the Green's functions will be 
shown to represent the free energy and the correlations of resolvent opera- 
tors resp. of the matrix model. The self-interaction of the scalar field and 
its interaction with background sources will be given as formal infinite power 
series, which receive 'short distance' (see below) corrections throughout the 
loop expansion of the effective theory. This latter expansion will be identified 
with the topological expansion of the matrix model. 
The unpublished work [TH] has been fused into the present article. 

The plan of the paper is as follows: In section 2 we collect some technical 
material concerning Hermitean 1-matrix models taken from the literature 
PP - [1], [LT]- (Experts may skip this section.) Section 3 is devoted to the 
construction of the effective Lagrange in tree graph approximation which is 
identified with the leading large N order for the matrix model. In section 4 
we discuss the non-leading terms of the matrix model topological expansion, 
i.e. the loop corrections in Lagrangean parlance. We end in section 5 with 
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some conclusions. 



2 Hermitean 1-matrix models 

The partition function of Hermitean 1-matrix models (pQ, [2]) is given by a 
(potentially formal) integral over Hermitean matrices 



Z N (t)= I e~ NtlV M dM , 
dM Yl dM a ] I Re dM ij } I Im dMi 



13 , 

i<j i<j 



Ma are real variables; the potential V(M) is given by V(M) = J2 n >i tnM n , 
with t n denoting some coupling constants. Ensuing averages are introduced 
as usual: 

(f(M)) = \\ f(M) e - NtTV{M) dM. (2) 



Z „ 

The free energy T related to the partition function by Z^it) = e r can 
formally be expanded in powers of -k (the so-called topological expansion) 

jr = Y^N 2 - 2h T {h) . (3) 

h=0 

Following well established paths [3] we introduce a derivative called the 'loop 
operator' 

— = -f— - (4) 

dV(z) Z3+ 1 dtj ' 1 ' 

z denoting here a complex number. Connected correlators of the resolvent 



N z-M 

are obtained by multiple application of the loop operator to the free energy 

w.{«, ...,«.) = jv— (w M ■ ■ ■ ffWU = WWW)- W^f- 

(5) 

The subscript 'conn' refers to the connectedness of the correlator. 
The correlation function W n inherits from the free energy a topological ex- 
pansion 

00 

W n ( Zl ,..., z n ) = £ N- 2h W^ h \ Zl , ...,z n ). (6) 

h=0 

2 



It turned out that the most efficient tool for the determination of the partition 
function (the free energy resp.) and of the correlators of resolvents are the 
so-called loop identities, that is, Ward identities for invariance of the defining 
integrals under reparametrisations of the matrix model integration variables 
of the type 

SM = e—^— 

z-M 

with e denoting an infinitesimal real parameter and z an arbitrary complex 
number. The invariance of the partition function under this variation leads 
to the identity 

(W 1 (z)r + j- 2 W 2 (z,z) = ^(tr^). (7) 

The l.h.s. of the latter equation derives from the variation of the measure in 
eq. (CD) under the variation ([6]) whereas the r.h.s. derives from the variation 
of the integrand exp(— iV tvV(M)). Applying to eq. (J7j) the loop operator— 
or equivalently, starting from a more involved variation substituting instead 
of eq. (jnj) 

5M = e l -— TT tr l — - 

z\ - M 11 Zi-M 
j>i J 

one can deduce further identities for multipoint functions. Eynard |1] uses 
those identities to establish for the correlation functions and finally also for 
the higher order corrections to the free energy a recursive procedure giving 
rise to a trivalent graphical representation (cf. [5] in the latter context). 
Invariance under the unitary group U(N) may be exploited to reduce all 
7V 2 -dimensional integrals above to iV-dimensional integrals over eigenvalues 
of the Hermitean matrices. The partition function, eq. ([T]), in particular 
becomes then [1] 

N 

z N (t) = I n^n^-^) 2 ri e ~ w(A!) ' ( 8 ) 

i=l i<j i 

where we have skipped — for convenience — on the r.h.s. of ([H]) the volume 
factor of the unitary group. We assume (as PQ-jl]) that the eigenvalues are, 
for large values of N, concentrated on a finite number of compact intervals 
on the real axis, say 

s 

(J A x with A x = [a 2 i-i, a 2 i], ai < a 2 < . . . < a 2s 
i=i 
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where the separate intervals [a 2 i-i,a 2 i] are supposed to be spread in the 
neighbourhood of diverse minima of the potential V(A). The spectral density 
Pat (A), concentrated on the above intervals, becomes a continuous function 
p(A) in the large N limit. The expectation value of the resolvent reads in 
terms of this distribution as 

Wl (p) = J d \0t (9) 
and the loop equation ((7|) becomes 

r^( w ) = 0^(p)) a + ( 10 ) 

Imp — io JS Z 
c 

The contour C of the integral on the l.h.s. is supposed to encircle in the 
complex plane the support of the spectral measure, but not the point p. The 
second term on the r.h.s. of ( TTU1) drops out in leading order of the large N 
expansion and one is left with the quadratic equation for the leading part 

J 2m p — oj J p — oj 

c 

= V'ip) Wi 0) (p) - f d J' {p) ~ V ' {uj) p{uj) = V'{p) W}%) + /(p), (11) 
J p-w 

where f(p) denotes a polynomial of a degree one less than the degree of 
is given as solution of the preceding quadratic equation, 

wi°\p) = \v\p) - Ivivipw+mp), (12) 

the sign of the square root term being dictated by the stipulated asymptotic 
behaviour W[ (p) ~ ^, p — > c>o. The expectation value of the resolvent is in 
this way related to a hyperelliptic Riemann surface given by the equation 

y 2 =(V'(X)) 2 + 4f(\). (13) 

One is free to choose / s.t. the polynomial on the r.h.s. has a finite even 
number, say 2s, of simple zeros in accordance with the above assumption 
on the support of the spectral density — besides an appropriate number of 
double zeros, s.t. y reads as 

y = M(X)y, (14) 
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2s 

f = Y[(\- ai ), (15) 

i=l 

where y now represents a reduced hyperelliptic surface [3] of the same genus 
g = s — 1 as y. The polynomial M(A) can be seen by inspection of the loop 
equation ( TTTT) to be given by a Cauchy contour integral (cf. [3]) 

where denotes a contour encircling all branch points Oj , % = 1, . . . , 2s and 
A. (A) is found similarly 



wf\\) 



dw V'(u) 



2f 2m\-u J \Uf =1 ^-a J ) 



e£i(a 



where C now encircles the branch points but not A. 

To fix the positions of the branch points one has to impose particular bound- 
ary conditions. Following [I], [T7] we choose those to be given by 'filling 
fractions' 

1 I dX m 



2 J 2m 

Ai 



which are the relative weights of the various spectral intervals. The insertion 
of the topological expansion (jBJ) into the loop equation (jTTJT) gives rise to 



h-X 



(ic - 2w['\p))w[ h \ P ) = w[ hl \p)w[ h - h \p) + wt 1] (p,p), (17) 



c 

The solution of ( 1TTI) under the condition 

dx 1 W 2 (x 1 ,x 2 ) = Vj 



(stipulating that the filling fractions are to be chosen independently of the 
form of the interaction potential V(t)) is found, [I], to be 

Wl%) = Res^f^^- (]T W[ h ~ m \x)w[ m \x) + Wt 1} (x,x) 



-f x^a-i dp y(x) 



\m=l 
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with dS(p, x) denoting a meromorphic one form with respect to the variable 
p and a multivalued meromorphic function in x on the reduced Riemann 
surface ( fT5l) . dS is characterized uniquely by the following properties: 

dS(x, x') = — h finite, 

x^x' x — x 

dS(x, x') =_ - + finite, 

x—fx' X X 

(x denoting the hyperelliptic involution of x) where the quoted singularities 
are the only ones on the reduced surface, and 



dS(x, x') = V j = 1, . . . , s — 1. 
The explicit expression for dS(x, x') is 



s-1 



dS(x,x') = Vpfl -J— -J^CA^Hx) dx (19) 

1 f dx 1 -rr 

with Cj(x') = ® — ; and a(x) = TT(x — ai). 

' 2m J x /rtxjx-x> 1 ; f \ K l> 



Lj(x) denote polynomials which are constitutive parts of the normalized 
holomorphic 1-forms of the reduced surface, 

iOj(x) = — dx 
^ix) 



<j) Uj{x) = 2iriSji for i,j£l,...,s-l. 



One needs for the resolution of the recursion, eq. ( 1181) . in explicit form the 
two-point function W^ip, q). This is found as solution of the loop equation 
for the two-point function — which emerges from another application of the 
loop operator to both sides of eq. (JTTJ) — in terms of the 'Bergmann kernel' 
B(p,q) (cf. H): 
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B is a symmmetric bidifferential on the reduced surface with a unique second 
order pole at the coincidence point p = q, B(p, q) = ( gzgp + finite) for p — > g 
and satisfying the normalization conditions 

y B{p,q) = 0, j = 1, . . . ,s - 1. 

The Bergmann kernel is related to the above introduced differential dS as 
follows: 

B M = ~(^ + dS(p iq )^dq. (21) 

The solution of the recursion relation (|T8l) and its generalization for multi- 
point correlation functions in terms of dS(p,q), B(p,q), y(q), i.e. in terms 
of intrinsic quantities of the hyperelliptic surfaces ffT41) and ffToT) is due to 
Eynard. For details of this formalism we refer to the original work [4] and 

3 The planar approximation in terms of tree 
graphs of an effective Lagrange 

The determination of the multipoint correlators of resolvents in leading order 
of i may either proceed by inspection of the multipoint generalisation of the 
loop equation ffTTj) — this approach has been followed by Eynard [4] — or may 
be achieved by repeated application of the loop operator to the above noted 
result, eq. ( 1201) . for the resolvent two-point correlation function. That is, one 
applies the loop operator either to the loop equation or to its solution. We 
will follow here the second route. 

The variation of the Bergmann kernel, eq. (1211) . due to an infinitesimal 
change of the branch points of the underlying hyperelliptic surface is given 
by one of Rauch's variational formulas [18] : 

SB(p,q) = ^B(p,[a j ])B{q,[a j ])6a j . (22) 

3 

B(p, [a,j]) governs the asymptotic behaviour of B(p, q) in q near af 
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In the following we will use 

B(p,q) 1 1 



B{p,q) 



dp dq 2 (p — q) 2 

and _ 

dz B(p,z) 

2m (z - (Lj) 1 / 2 ' 



B(p, [ aj \) = 2j 



Rauch's formula is to be combined with an expression for s ^ p ^ ■ The latter 
is found by first noting the relation, [5] , 

and to match this with the asymptotic behaviour y(q) = y([ai\)y/q — di + 
0((q — aj) 3 / 2 ), which implies 



Sy(q) 



1 ,r 1 5dj 



SV(p) 

Eqs. (E3D and (El) lead us to 

<faj 2£(p, [a,]) 

sv(p) i/(N) 

which with eq. (1221) gives rise to 

2s 



(25) 



r(o)/ \ 5_ (o) , . _ Bjpx, [oj]) B{p 2 , [ai])B(p 3 , [a,]) 



^W*) = ^<W) = E y(N) 

We will use the notations 



(26) 



Z? / " | <A ^f- / ^ (27b) 

f 2mJ 2m (z - a^f+^z - a^+V 2 { ' 



ai aw 



, dz y(z) x 

= ^ (-^ (27c » 



r(0)/ 



1 We remind in this context the identities y(q) = —2Wi (q) + V'(q) and 
gy( p ) V'Cg) = — ( g _ 1 p )2 ; the latter being a direct consequence of eq. JD- 



Pi 




P2 P3 

Figure 1: Three point function W^°\ The vertex factor is 



It will turn out that the collection of quantities (27a)- (27c) is complete in the 
sense that the resolvent correlators and the free energy can be expressed — as 
will be shown — to all orders in the topological expansion by those terms. One 
should note that (27b) is well defined also for i = j, since the order of the 
contours does not matter due to the fact that B(p, q) has as single singularity 
a second order pole with residuum one. Straightforward calculations using 
eqs. ( 121)]) and (|2U|) lead to the relations 

^ = (Ai(q) + A 2 (q))(yf ji ) (28) 

(A 1 (?))(y/ |i ) = (2/ + l)^5°(g) (EE 

{A 2 (q))(y f>i ) = -B{(q) 



5V(q) 



(A 3 (q) + A 4 (q))(B{(p)) 
(A 3 (q))(B{(p)) = (2f + l 



2s r/.°rO/ 



(30) 



6B f ' 9 

^ = (A 5 (q) + A 6 (q)+A 7 (q))(Blf) 

(A 5 (q))(Bif) = (2/+i)-y — ^ dm 

B f ' 9+1 B°( a) 

(A e (q))(Blf) = (2g + 1) ^ ~ » 
(A 7 (?))(S{f) = £ fcW . » 

The application of an operator Aj(g), j = 1, . . . , 7 to expressions not noted 



(29) 
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above gives zero. It can also easily be verified that 



and 



are 



SV(q) 6V(p) 

for arbitrary arguments p and q commuting operators acting on the basis 
(27a-c) (as it should be indeed). The 4-point function is found by applica- 
tion of (cf. eqs. (EE} and @2}) to wf. One finds: 



Wi 0) (pi, 



= £ (tffeOBffe) (£) iff (^-) 



,P4) 



SV( Pi ) 



+ S?(Pi)S?(p 4 ) 



1 

yi,< 



5 



0.0 



1 

4 2s 



B?fei)B°fe) 



i=l \ fl,*/ r =l r =l i=l \ y l>«/ t=l 

Eq. (|2T)]) and flHTl) designate the starting for an effective Lagrangean descrip- 
tion of the resolvent operator correlation functions. The most obvious hint 




+ (p 2 < > p 3 ) +(p 2 < > p 4 ) 



+ (p 2 < > p 2 ) +(Pi< > p 3 ) +(Pi< > p 4 ) 



Figure 2: Various contributions to the 4-point function, (a) Second order 
contributions of the cubic interaction. (b) + (c) Genuine quartic interactions 
without and with a derivative field, the latter being indicated by a bar. The 
various vertex factors are not accounted for. 



from both equations is that the interactions are concentrated at the branch 
points a, , i = 1, . . . , 2s and that an interaction point is connected to the 
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external resolvent operator by a 'propag ator' 5?(p). 3 The 3-point correlator 
can obviously be related to a cubic interaction: 

Sf = f 4°>M d x with £?>(*) = -1*®. = l£ (32 ) 



j/i(ar) 3! j/i 3! 
with the correspondences 

[dx = jr\ , £(^1(p)=b°(p). 

I \x=CLi 

J 1 = 1 

The first bracket on the r.h.s. of (l3"Tj) represents then the second order contri- 
bution of the cubic interaction to the 4-point correlator where one propagator 
is exchanged between the interaction vertices at places dj and aj. The second 
and third term on the r.h.s. of (1311) are generated by quartic interactions as 

(o)(«) _ 3 VA V A _ ^ ( %2,A ¥>V 



and 



^-/^(-fj^-E^ffJ^ (33) 



S? m = [ dx'^ = ± 1^1^ (34) 
J yl{x) 3! 1! ^ y 2 lti 3! 1! 

where the derivative (dcp)i = (dcp)(ai) gives rise to the factor Bj(p ). The 
Lagrange resulting from these quartic interactions is 



ViJ 4! \y-{] 3!1! 

We will have to make use of couplings with derivatives of arbitrary order: 
B{ (p) will point to an f-th order derivative field at place a,, 



(3fy)< <p(p) = S/(p) (df ( p) i (d'(p) j = B(f. 

Consecutive application of loop operators to W± (pi,P2,P3,Pi) will generate 
the (connected) correlators of an increasing number of points. The applica- 
tion of the pieces A 4 fl29b) and A 7 (!30b ) lead to the creation of new vertices 
inserted in external and internal propagator lines resp., whereas the applica- 
tion of the remaining parts of gives rise to a diversification of the already 



2 We will treat the Lagrangean formalism on a formal level without making concrete use 
of the fact that the Bergmann kernel B(p, q) is the Green's function of a chiral derivative 
field on the reduced hyperelliptic surface as it was noted and used by Dijkgraaf and Vafa 
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existing vertices. Let us concentrate for a moment on those contributions to 
the n-point function emerging from a single vertex, that is the part of the 
amplitude which will be related to the n-field terms of the effective Lagrange. 
There will be terms with 3 + k, k — 0, . . . , n — 3 propagators without deriva- 
tives out of the total number n emerging from an interaction point — and 
there will be corresponding pieces in the Lagrange with the same number 
of fields without derivatives. Keeping in mind that the contribution to the 
n-point function going along with the n-field vertices has to be for itself sym- 
metric in the arguments pi, . . . ,p n one may proceed as follows: One first of 
all generates k propagators without derivatives by applying k times Ai to 
the 3-point function, eq. ff26l) . and afterwards A 2 for another n — (3 + k) 
times to generate the propagators with derivatives. The symmetrisation in 
Pi, . . . ,p n has to be performed afterwards. That is, the 1-vertex contribution 
to the leading order n-point function is given by the symmetrisation of the 
following expression 

n— 3 n— 3— k k 

I! A2(P3 +fc+ r)n A l(P3+t) WP( PuP2 ,p 3 ). 
k=0 r=l t=l 

One easily extracts therefrom the complete effective Lagrange beyond the 
two lowest orders (eqs. (J3"2j) - (}3"1|) ) 

oo 2s 2s 

S (0) = J> = ££4 0) (a,) - (35) 

n=3 n>3 i=l i=l 

as a formal infinite series. For the sake of illustration we quote the concrete 
expressions for C$ and 

5 V Vl Vi) 5! + l yf) 4!1! + {yf) 3!2! + { 6 yf) 4!l! 

(36) 

.yl , ,, n 2/22/3 i nc yA v? 6 , do<-yl my-A 9^ 



4 0) = f-4054 + 450^ - 105^ ^ + ^4 - 60^| 

6 V yl vl vl) 6! V yl yl 



5!1! 



+ , y — iiiiii— + l 15 5?J-5iir- (37) 
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We state our first main result as 
Theorem 1: 

The connected correlation functions of the resolvent operators are in leading 
order of -k given by the tree graphs of the Gell-Mann-Low series of the effec- 
tive Lagrange C^°\ eq. (1351) . For the determination of an n-point correlator 
one has to evaluate up to the n-th order, £( 0,ra ) = ^fc=3^i°^' an d to 
specify all n-point tree graphs deducible from £(°' n \ 



The proof of the theorem is a proof by induction. Suppose that the state- 
ment of the theorem has been found to be true for correlators with up to n 
points. Applying to the n-point function a further time the loop operator 
one generates either new triple vertices inserted in all possible internal and 
external lines by the action of A7 (eq. (130b )) and A 4 (eq. (|29b )) or one adds 
a new external line to any of the already existing vertices by the action of 
A 1; A 2 , A 3 , A 5 and A 6 . 

Some thought reveals that in this way all connected tree graphs with (n + 1) 
end points related to = Y^^l Cf^ are generated if one takes the 

induction assumption for grantedjj 



4 Loop corrections 



4.1 1 loop order 

Starting point is the 1-loop version of eq. (|TH|) : 

W<%) = £ Res J_ wf (*,*). (38) 

' ' x—*n- an i/lrl 



+(h dp y(x) 



i=i 

dS(p, x) and its derivatives can be expressed at the branch points x = a r by 
Bf(p) as follows: We define 

dS(p,[x]i) = 2dS(p,x^ 



y/X — 0>i 

B(p,[x]i) = 2 B(p,x)y/x - a~i 

y([x]i) = y(x) 

yX — (ii 



3 One should note in this context that the factorials of the Gell-Mann-Low series are 
completely absorbed, as there are no symmetry factors left in tree graph order. 
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and recalling eq. f[2Tj) one is immediately lead to 

~ ldS(p,[x]i) 1 dS(p,[x]i) , , 

B(pMi) = l 1 V + 7>(x-a i )d x — . (39) 

Differentiation of the last equation with respect to x and then putting x = 
gives 

x lx=a * dp 2k + 1 1 v ' v ' 

From the inspection of the concrete appearance of W2(p,p) — that is, the form 
of the Bergmann kernel, of eqs. (fl""9l . f[2T]) resp. — we obtain 



W (0 \x,x) = if - if J — + -Y 

16 ' (x — di) 2 16 ' x — a* x — a,- 4 ^— ' x — a; 

ly 1 ly 1 1 | ly PjK) 

16 ' (x — ai) 2 8 ' aj — a~ x — a; 4 ^— ' x — a.; 

j=l V * 7 i,J=l 1 3 % i=l 4 

2s ^ 2s D 0,0 

= -y (4i) 

16 ^ (x - ai) 2 4 ^ x-di K J 

i=l v ' i=l 

where the asymptotic relation (p, p) \ p ^oo ~ 1/p 2 is used for the second 
equality and the last equality can be inferred from a straighforward evaluation 
of the double integral constituting B®f. The polynomials P m are related to 



the polynomials Lf 



dz 



3=1 



Pm{x) = - y^Lj(x) <b -j= — m = l, ...,2s. 



A, 



(z - a m )y/cr(z) 



Inserting ( HIT) into ( 1381) and taking into account ( l40~i) we arrive at 

The first two terms on the r.h.s. of ( [4*21) emerge from the double pole part 
in eq. (j4"Tl) . whereas the third term is due to the part in (j4"Tl) with a single 
pole. The latter contribution represents the piece which could have been 
anticipated — taking the Lagrangean point of view seriously — as a tadpole 
correction of the cubic interaction = to the 1-point function which 
is graphically depicted in fig. E^c). The repeated action of the loop operator 
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o 



p p p 

(a) (b) (c) 

Figure 3: Graphical representation of the 3 terms contributing to W[ . The 
small circles represent vertex insertions with topological index 1. 



Figure 4: 1-loop corrections to the 2- and 3-point correlation functions 



on the tadpole contribution gives rise to the complete set of 1-loop graphs of 
£(o) = ^2 n>3 Cn^ (eq. f[3"5l) ). some of those being depicted in fig. HI 
The first two terms on the r.h.s. of eq. (|42|) may be viewed as 'short distance 
corrections' to the above tadpole contribution. To characterize those terms 
and their generalisation to be introduced below we ascribe (in a slightly ad 
hoc fashion) to the quantities Bf , and yi+x',i the mass dimensions x+1/2 and 
x' resp., which amounts for the former quantity to the attachment of a mass 
dimension x + 1/2 to a derivative field d x ip. We also introduce the notion 
of a topological index h(V n ) for an interaction vertex V n of n emanating 
propagators 

h{V n ) = ±deg{V n ) -~n + l (43) 

(and the analogous notion for the corresponding part of the interaction La- 
grange) where the mass dimension deg(V^) counts the sum of indices related 
to ends of propagators emanating from V n and the sum of dimensions of fac- 
tors y, t . attached to V n . The topological index of vertices of the tree graphs 
considered in the preceding section and therefore also the Langrange 

£(0) 

(eq. ( 1351) ) is according to this definition zero. The same holds for the ver- 
tices of the 1-loop graphs displayed in figsE^c) and HI The two short distance 
corrections in (1421) on the other hand are of topological index 1, as well as 
the first non-leading correction to the tree graph Lagrange which is obtained 
by application of Ai, A 2 , A 3 , A 5 and A 6 (cf. eqs. ( jZEHSDD ) to those terms in 
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question. One finds in this way in particular 

00 

£« = (44) 



(i) _ /3y| 5y 3 \ <^ 2 / 1 j/2 \ (g (gy) 



=i 



1 2/2 \ v> , ( i i \ c^v 2 



^ 2 V4yf g^J 2! + V 4^y 111! 

/ 1 i \ (g y ) 2 Ai u (ay) 

^24y^ 2! ^VSl/iV 

-(l) / yf 105^3 35t/Ay? 3 ( y\ 15y 3 \ ^ 2 ^ 

3 " V Vl 8 3! V yf 8^ 2!l! 

3y 2 \<p(d<pY , / 3y 2 W(«9V) , /l 1\(<^) 3 



' AyfJ 1!2! ' V 2!1! Vl2y 3 

1 J_\ y (fly) (,9V) + /5 1 ^ y, 2 (9V) 



v 4^y mil! \SyfJ 2!1! 

To summarize the preceding: The next to leading order contributions to the 
connected correlation functions are given by expressions corresponding to 
1-loop graphs with vertices of topological index zero or to tree graphs with 
one vertex from eq. (j44l with topological index 1. 



4.2 More loops 

4.2.1 A preparatory calculation 

Proceeding from order (ft, — 1) to ft in the topological expansion one has to 
resort to the loop equation ( fT8l) . An immediate consequence of the recursive 
nature of equation ( 1T81) is that the 1- and 2-point functions have the general 
form 

W ( ?\x)=Y J B f i {x)u { tXiJ) (45) 

and 

W^(x,x) =^B{{x)Bj{x)4 h \i,f-,j,g) (46) 

id 

f,g 

where the functions uj^ do not depend on x. Confronting eq. (ITBl) with 
the last two equations we see that we have to deal with the evaluation of 
expressions 

^^x) 1 , 

dp y(x) 3 k 
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We do this for arbitrary integer values / and g and the special (most involved) 
case i = j = k, merely quoting at the end the results for the other (simpler) 
cases, which are practically contained in the case of coincinding arguments. 
Cauchy's integral representation for the residuum and the representation eq. 
(27a) for B[ and Bf lead to 

Res *S^*B!{x)Bl(x) 
x^a, dp y(x) 

dx dS(p,x) n / dw B(w,x) „ / dv B(v,x) 



2ttz dp y(x) J 2m (w - a;)/+V 2 J 2m (v - a^+V 

The contour in x is by definition positioned outside the contours of those for 
the variables w and v. The order of the latter contours is immaterial but we 
assume for the sake of definiteness that the contour in w surrounds that in 
v . Pushing the x-contours in (I47p through the other two one picks up two 
additional terms from the double poles of B(w, x) and B(v , x) at w = x and 
v = x resp., 

(H7D = 4 I— - I— - I^ dS ^ x ^ B(w x)B(v x 

&i di &i 

, dw 1 / dv 1 /l / dS(p,w) ~ 

P 2m(w- ai )f+^ J 2m(v-a t )9+V 2 \2 w \dpy(w) {W,V ' 



ai ai 

+ 



1 f dSj^v) - \\ 
2 d \-dp-y^ B ^)) ' 



The first term on the r.h.s., let us denote it as ff48a ). where the x-contour is 
now the most inner one, can easily be computed since the integrand has a 
simple pole at x = ai. Taking into account eq. ( )40l) we obtain 



lBMB[fB?f 



2 yi,i 

The two other terms on the r.h.s. of (HHl) . denoted (TJEb) and (HHb ) resp., 
become after partial integration 



Eb) = (2/ + 1) 



dw 1 I dv 1 dS(p, w) B(w, v) 



2ni (w — ai)f +3 / 2 J 2m' [v — a^^ 1 / 2 dp y{w) 

ai 

(49) 
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and 

dw 1 / dv 1 dS(p,v) B(v,w) 



2m (w - ai)f +1 / 2 J 2m (v - a^+^l 2 dp y(v) 



a, 



(50) 

The f-integration in the last equation — being now the most inner one — can 
immediately be executed and one finds with eq. ( HOI 



dSOD = (g + \) E WZi^ B v z *>* withz ^ ~ d '~ 



t\ y([x]i 

k+ r +t=g+l yu 1 

(51) 

The evaluation of proceeds by pushing the w contour through the v- 
contour. One obtains in this way one term analogously to eq. (151jl and an 
additional term from the pole at w — v yielding: 

m = u+\) E w^ B}(p ') B S z *>* 

k+r+t=f+l 

+i(2/ + l)(2p+l) E ^-l B i(p)Zi,i. (52) 

k+l=f+g+2 

Putting pieces together, i.e. eqs. (HHk). ( 15TT) and ( 152]) and taking care of the 
cases of non-coincident indices, one arrives at 

Res ^pP2A B Ux)Bl{x) 
x^ ai dp y(x) 3 k 

+ Sk,i(9 + \) E 2^TT s * r (p)^f^ 

r+m+t=g+l 
r+m+t=f+l 

+ 5^6^(29 + 1)(2/ + 1) £ _L_5 [(p )z M . (53) 

r+Z=/+</+2 

It is worth noting for later reference that from all terms on the r.h.s. of eq. 
( 1531) only the last term is responsible for the creation of a vertex with an 
increased topological index (The index of the fused vertex is then equal to 
the sum of the two (non-fused) original vertices if it results from the first 
term on the r.h.s. of eq. (118j) or it increases by one unit if it originates from 
the second term.) 
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4.2.2 Free energy at higher loop orders 

We adapt momentarily our notations to those of ref. [5] by writing Z(t) = 
J dMexp(-gtrF(M)) = e T with V{M) = J2 n >i l rM n and remind that the 
occupation numbers were above designated by e». Chekhov and Eynard [5] 
noted that the scaling relation for the free energy 

with h = j^, is synonymous with a formula which gives as 'integral' of 

W^. The details of the integration operator, called H in [5], is immaterial 
for our purposes. We have only to note that H applied to B gives 

B(; P ) 1 

^■^ = -2^ 
to find with the help of fT4"5| and the previous scaling relation 

F (h) = -2fh2 H *' W ™ {9) = 27^2^^ )(M) ^ (54) 

(>i 

where the I = part of eq. (14"5| drops out since it leads to a vanishing 
residuum. Eq. (j54p holds for h > 2. What concerns the expressions for 
and J-(°' we refer to [H] and [H] and in particular [21] . 

4.2.3 Short distance corrections 

We want to work out the repercussions of the last term in (1531) on the higher 
order corrections to the free energy, the 1-point function and the effective 
Lagrange resp.. We restrict here our attention to 'local' contributions, i. e. 
to those quantities which only depend on functions y x ^ and Bf from one 
and the same branch point at a time and do in particular not depend on 
propagators Bff, Bfj y connecting the same or different branch points. We 
use here and in the following a hat to designate 'local' quantities. Starting 
from the short distance corrections to W[ (jp), the two first terms on the 
r.h.s. of eq. (EE]) . 
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one obtains a 'local' contribution to W% (p, q) by applying the parts Ai, A 2 
and A 3 of ^ to W^: 



i=i 

2s 



Inserting the local part of eq. (j4"2"|) and (1531) into the loop eq. (TT5T) and 
restricting to the last part of eq. (1351) one finds 

A/ 21 ^ 29^3 35 y A 
+ ^l v -32^ + 32^T-128^J t ^ (Pl) 

^ / 63 29 y 3 \ x 



i=l 
2s 



and from eq. ( 1561) a local contribution to the free energy 

' 21 yl 29 y 2 y 3 35 y 4 
160 yf 128 384 



^ ,2) = £ ( 

i=i ^ 



(57) 



Applying the loop operator to W[ (p) one finds W\ [pitPvi ^3 

etc. and can extract from this the two loop short distance correction to the 
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effective Lagrange 



4 2) = 





1323i/ 2 5 495y|/3 _ 2205y|/4 _ 8175^1 63^5 1785 y 3 y 4 
" 32 y? 4 yf 32 y x 7 128 y 7 2 yf 64 yf 
^ ' '441 225 y 2 2 y 3 385 y 2 y 4 725 yj 10b y 5 \ <p (dtp) 
K 32yf 8 y 7 32 yf 128 y? 32 yfj 1!1! 

63 yf 145 y 2 y 3 
16 y\ 32 yf 
g_(gV) /63y 2 2 29y 3 s 
1!1! \S2yf 32 

W (gV) , f 

-yi/ HI! V 80y5 y 

(gg) (gV) / 29 1 \ (gV) (gV) 
1!1! + V 128y^ 1!1! 

It should by now be clear how to proceed to the determination of higher 
order local quantities. Let us assume that W[ (p), . . . ,W[ (p) have been 
determined. One obtains (p, q) by loop differentiation of and then 
W[ h+1 \p) by insertion of YL=i w[ h+1 ~ m \x)W[ m \x) + w!f\x,x) into the 
last term of eq. ( 1531) and therefrom !F^ h+1 \ One arrives in this manner for 
example at 

~ (3) = /2205 y\ _ 8685y|y3 15375yjy| 5565y|y4 _ 72875 yg 

^ V 256 vl° 256 y? 512 yf 256 2/1 8 2 i504y 1 7 

5605 y 2 y 3 y 4 _ 3213 yfys 21245 yf 2515^ 5929^6 _ 5005 yr\ , . 
256 y 7 ' 256 y 7 9216 yf 512 yf 1024 yf 3072 yf// 1 J 

In the table of all single Lagrange functions two pieces in h = are 
missing because we start with W 7 ^ ^. We absorb from 1-loop order onwards 
the function W± as into the Lagrange and get finally: 

£(°) = 4 0) + 4 0) + ... 

C (h) = c W + c W + c W + c (h) + ^_ fo r /i>l. 

One may also introduce h > 2, as the local part of the free energy. 
What concerns the full free energy Fl h \ cf. eq. (I54j). 



4 W r -[ ' ) and instead supply the ingredients of the above introduced graphs, the 

vertex factors and the propagators. 
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To take the two special cases appearing in the Lagrange also out of the 
correlators, one has to define 



W k = W k -8 ktl W™ -6 kt2 W™. 

The expression (0| . . . |0); i 00 ps, conn denotes the sum of all I loop diagrams of 
(0| . . . |0)„ 

The main result of the present paper is comprised in 
Theorem 2: 

With the Lagrange £ = C^°' + jjp £^ + -m O- 2 ' + • • • the Hermitean 1-matrix 
model correlation functions are given by 

00 ^ 

W k (p!, . . . ,p k ) = J2j^l < I ' e E -i £(fll) |0), loops, conn- 

1=0 

For the evaluation of the Gell-Mann-Low series the propagators 

I I I I 

(df(p)i (p{p) and id^<p)i id 9 (p)j have to be taken as B[ (p) and B- f resp.. 

sketch of proof: One has first of all to show that the interaction vertices 

are all generated by the repeated action of the loop operator on w[ h \ This 
can be proved inductively: Assuming that it is found true to (n — l)-th order 
the step to n-th order is done through the loop equation ( fl8l) . It is in this 
context particularly convenient to consider the free energy F (h) . All terms 
of eq. ( |53i) besides the last term on the r.h.s. of eq. (1531) (and leading to the 
local correction J r ^) have a graphical appearance as shown in fig. [5j Those 




Figure 5: Graphical appearance of terms contributing to the free energy 
not being of the local type T^> 

contributions to the free energy depicted in fig. [5J are related to contribu- 
tions to ip) of fig. E] via action of loop differentiation. The amplitudes 
shown in fig. E] are uniquely connected by the loop equations to lower order 
terms given in fig. [71 (We select for the purpose of our argument a term 
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Figure 6: Contributions to W ] 




Figure 7: Lower order amplitudes related to those of fig. E] 



with k ^ i,j which does not give rise to short distance corrections at the 
extra vertex placed at a*,.) But for the latter amplitudes holds the inductive 
assumption and therefore also for all contributions to J-^ which are not be- 
longing to the local part of J r ^ h \ It is then in particular true that the 
local vertices making up are those which were already present in JF^ -1 ) 
(covered by the inductive assumption) and the new local vertex for J 7 ^' (and 
then also for W^^C^i etc.). This concludes the argument. 
The missing part for a full proof of the above theorem, that is, a control of 
symmetry numbers in 2 loop ordeiEl and beyond, requires a detailed com- 
binatorical analysis and will be given in a separate forthcoming publication 

5 Summary and conclusions 

The purpose of the paper was to develop a Lagrangean formalism of Her- 
mitean 1-matrix models taking as starting point Eynard's technique [4] en- 
coded in the hyperelliptic Riemann surface associated as spectral surface to 
the matrix model. Eynard's approach is based on 2 ingredients: 1.) The dif- 
ferential 1-form attached to the defining equation of the hyperelliptic surface 
and 

2.) the two-differential, the Bergmann kernel, on the reduced surface. Those 
ingredients are found to show up in the solution of the loop equations for the 
connected correlation function of resolvent operators given by an iterative 

5 The proper appearance of the 1-loop order has been checked. 
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nested system of Cauchy contour integrals. We resolved this nested system 
into the rules for an effective Lagrange of a scalar field propagating on the 
reduced surface, s.t. the propagation is by the Bergmann kernel with multi- 
ple self-interactions of the scalar field taking place at the branch points. The 
interactions are represented as a formal infinite power series in powers of the 
scalar fields and its derivatives and are subject to short distance corrections- 
due to the singularity of the Bergmann kernel at coinciding points — at all 
orders of the topological expansion. Our motivation for the endeavour to con- 
struct a Lagrange formalism is to find a new approach to critical behaviour 
via the renormalization group. The test of the usefulness of the formalism in 
this context is still ahead of us. 

Evidently, the formalism reveals an abundantly rich true Lagrangean struc- 
ture behind Eynard's trivalent graphical representation. 
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